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Abstract.  Numerical  and  experimental  investigations  are  performed  for  the  rarefied  gas  flows  in  pumping  channels  of  a 
helical-type  drag  pump.  Modem  turbomolecular  pumps  include  a  drag  stage  in  the  discharge  side,  operating  roughly  in 
the  pressure  range  of  lO'^-lOTorr.  The  flow  occurring  in  the  pumping  channel  develops  from  the  molecular  transition  to 
slip  flow  traveling  downstream.  Two  different  numerical  methods  are  used  in  this  analysis:  the  first  one  is  a  continuum 
approach  in  solving  the  Navier-Stokes  equations  with  slip  boundary  conditions,  and  the  second  one  is  a  stochastic  particle 
approach  through  the  use  of  the  direct  simulation  Monte  Carlo  (DSMC)  method.  The  flows  in  pumping  channels  are 
three-dimensional  (3D),  and  the  main  difficulty  in  modeling  a  3D  case  comes  from  the  rotating  frame  of  reference.  Thus, 
trajectories  of  particles  are  no  longer  straight  lines.  In  the  present  DSMC  method,  trajectories  of  particles  are  calculated 
by  integrating  a  system  of  differential  equations  including  the  Coriolis  and  centrifugal  forces.  Our  study  is  the  first 
instance  to  analyze  the  rarefied  gas  flows  in  rotating  frame  in  the  presence  of  noninertial  effects. 


INTRODUCTION 

Today,  molecular  drag  stages  are  widely  used  in  conjunction  with  a  turbomolecular  pump  in  order  to  improve  the 
pumping  performance  of  it.[l]  For  example,  compound  or  hybrid  molecular  pumps  have  gained  wide  attention 
because  of  their  increased  pumping  throughput  and  discharge  pressure.  The  drag  pump  may  be  of  the  Holweck 
(helical-type)  or  Siegbahn  (disk-type)  design.  In  the  present  study,  the  pumping  performance  of  a  helical-type  drag 
pump  (HTDP)  is  numerically  and  experimentally  analyzed.  The  flow  within  a  pumping  channel  of  a  HTDP  varies 
from  the  molecular  transition  to  slip  flow  regime.  Two  different  numerical  methods  are  used:  one  is  a  continuum 
approach  in  solving  the  Navier-Stokes  (N-S)  equations  with  slip  boundary  conditions,  and  the  other  is  a  stochastic 
particle  approach  through  the  use  of  the  direct  simulation  Monte  Carlo  (DSMC)  method. 

Panos  et  a/.  [2]  cfefined  the  separation  pressure  between  the  viscous  and  the  free  molecular  state,  and  they 
extended  Gaede's  equation,  which  was  used  for  the  helical  single-groove  drag  pump  in  the  free  molecular  flow,  to 
the  helical  multi-groove  drag  pump  in  both  the  viscous  and  free  molecular  states.  Nanbu  et  a/.[3]  numerically 
studied  the  moleculartransition  and  slip  flows  in  helical  channels  by  using  the  DSMC  method.  Sawada[4]  analyzed 
the  pumping  mechanism  by  solving  the  N-S  equations  with  slip  boundary  conditions.  However,  the  previous 
studies  [2-4]  were  limited  to  the  quasi  three-dimensional  (Q3D)  model  of  Couette  flow.  Namely,  the  three- 
dimensional  (3D)  rotating  helical  channel  was  considered  as  a  rectangular  groove  facing  a  wall,  which  is  linearly 
moving  along  the  groove;  see  Fig.  1 . 

In  the  recent  study  of  Niu[5]  and  Cheng  et  a/.  [6],  the  rotating  helical  channels  were  three-dimensionally  modeled 
to  investigate  the  curvature  effect  of  channels,  and  a  laminar  compressible  N-S  solver  with  first-order  velocity  slip 
condition  was  used. 

In  the  present  study,  the  DSMC  method  using  the  3D  model  is  employed  to  investigate  the  effect  of  Coriolis  and 
centrifugal  forces  on  the  pumping  performance.  Also,  the  formulation  of  second-order  slip  boundary  conditions 
(SBC)  is  utilized  in  the  N-S  method.  As  a  validation  for  our  work,  the  present  numerical  results  for  Q3D  and  3D 
models  are  first  compared  with  the  previously  known  numerical  and  experimental  ones  obtained  by  Nanbu  et  a/. [3] 
Of  the  two  models,  the  3D  model  gives  more  accurate  results  than  does  the  Q3D  model.  Our  experimental  results  are 
presented  in  the  pressure  range  of 0.02- 4  Torr  and  are  compared  with  numerical  ones  for  both  models.  Also,  we 
assess  the  validity  and  usefulness  of  the  N-S  simulations  in  the  moleculartransition  regime. 
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COMPUTATIONAL  MODELS 


Two  different  kinds  of  helical  rotors  are  Illustrated  in  Figs.  1(a)  and  1(b),  where  a  is  the  helix  angle,  a  the 
channel  width,  b  the  ridge  width,  d  the  channel  depth,  and  At/  the  clearance  between  rotor  and  stator.  The 
Knudsen  number  is  defined  as  Kn  =  XI  d  ,  where  A  is  the  mean  free  path  of  molecules.  The  geometrical  dimensions 
of  channels  are  given  in  Table  1.  The  channel  depth  of  rotor  1  is  4  mm,  while  that  of  rotor  2  is  continuously  varied 
from  8  mm  (at  inlet)  to  4  mm  (at  outlet).  Also,  the  channel  shape  of  rotor  2  consists  of  two  parts,  i.e.,  a,  =  40°  for 

the  upper  part  of  channels  and  a2  =  15°  for  the  lower  part  of  channels,  respectively. 

Two  different  numerical  methods  are  used  in  the  present  investigation.  The  first  one  is  the  continuum  approach, 
in  which  the  N-S  equations  with  slip  boundary  conditions  are  solved  numerically.  The  second  one  is  the  molecular 
approach,  in  which  a  direct  physical  simulation  is  performed  by  considering  interactive  motions  of  molecules.  The 
DSMC  method  of  Bird[7]  is  employed  to  solve  the  problem  of  internal  rarefied  flows  in  pumping  channels.  To  save 
space,  we  do  not  discuss  here  the  procedures  in  obtaining  solutions  of  the  N-S  equations;  but  we  refer  the  reader  to 
Sawada.[4] 


DSMC  SIMULATIONS 

The  simulated  physical  space  is  divided  into  a  network  of  cells .  Figures  2(a)  and  2(b)  illustrate  the  computational 
grid  of  the  3D  channel  for  rotors  1  and  2,  respectively.  A  rectangular  grid  system  is  used  in  Q3D  model,  as  seen  in 
Fig.  3.  In  the  case  of  3D  model,  each  cell  is  subdivided  into  five  tetrahedral  subcells,  which  have  four  well-defined 
triangular  faces,  to  enable  tracking  the  molecules  from  one  cell  to  another. 

A  large  number  of  molecules  is  uniformly  distributed  in  the  computational  domain,  and  the  initial  velocities  of 
molecules  are  evaluated  from  the  Maxwellian  distribution  function.  The  gas  flows  entering  the  open  boundary  of  the 
computational  domain  are  initially  assumed  to  be  in  equilibrium,  with  zero  bulk  velocity.  The  bulk  velocity  is 
updated  by  iteration.  If  a  molecule  collides  with  the  wall,  a  new  reflected  velocity  is  assigned  by  the  surface 
boundary  condition.  Both  wall  and  gas  molecules  passing  through  the  pumping  channels  are  assumed  to  be  in  an 
isothermal  (300  K)  state.  In  the  present  simulation,  the  gas  is  nitrogen.  The  null-collision  technique  of  Koura[8]  is 
adopted  to  determine  the  number  of  collision  pairs  selected  within  a  cell  during  time  step  A/ .  The  collision  process 
is  accounted  by  using  the  variable  soft  sphere  molecular  model  of  Koura  and  Matsumoto.[9]  The  flow  field  is 
sampled  every  5  time  steps  during  20000  time  steps  after  reaching  the  steady  state. 

TABLE  1 .  Geometrical  dimensions  of  helical-type  rotors 


Rotor 

1 

2 

Rotor  diameter 

D 

168.6  mm 

168.6  mm 

Rotor  axial 
length 

L 

132.0  mm 

Lj  =  66.0  mm 

Z-2  =  66.0  mm 

Helix  angle 

15" 

a,  =  40" 

a 

o 

II 

Channel  width 

13.6  mm 

a,  =  42.2  mm 

a 

a^—  13.6  mm 

Ridge  width 

5.4  mm 

b^  =  5.4  mm 

b 

b^  =  5.4  mm 

Channel  depth 

d 

4.0  mm 

=  8.0  mm 

d^  =  4.0  mm 

Clearance 

Ad 

0.9  mm 

0.9  mm 

No.  of  grooves 

N, 

7 

7 

Rotational  speed 

rpm 

24000 

24000 
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(a)  (b) 

FIGURE  1.  Geometry  of  helical-type  drag  pump  (Q3D):  (a)  rotor  1;  (b)  rotor  2. 


(a) 


(b) 


FIGURE  2.  Computational  grid  system  for  3D  helical  channel:  (a)  rotor  1;  (b)  rotor  2. 


FIGURE  3.  Computational  grid  system  for  Q3D  model  (rotor  1). 
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The  main  difficulty  in  modeling  a  3D  case  comes  from  the  rotating  frame  of  reference.  In  particular,  trajectories 
of  particles  are  not  straight  lines  in  the  rotating  frame.  Thus,  trajectories  are  calculated  by  integrating  a  system  of 
differential  equations,  including  the  Coriolis  and  centrifugal  forces.  The  fourth-order  Runge-Kutta  algorithm  is  used 
to  integrate  the  following  equation: 


Ur  =  ,  (1) 

in  which  is  the  relative  acceleration,  and  is  the  trajectory  and  the  relative  velocity  of  molecules,  and 
(O  is  the  angular  velocity  of  the  rotor. 

From  Eq.  (1)  we  can  find  positions  and  velocities  of  particles  at  any  time  on  the  y-z  plane;  see  Fig.  2(a).  The 
component  of  the  particle  velocity  along  the  x-axis,  c^,  is  independent  of  the  frame  of  reference,  i.e.,  is  the 
same  for  both  moving  and  nonmoving  frames  of  reference. 

At  the  end  of  the  simulation,  the  dimensionless  pumping  efficiency  w ,  pumping  speed  S  {Us),  and  throughput 
Q  (Torr- //s  )  are  calculated  by 


w  = 


(2) 


S  =  w  A^  -^JSRT In  Kis^)/ 4.0, 


(3) 


where 


Ar(5,)  =  exp(-5|^)  +  JnSj  [l  .0+erf  (5, )] ,  5,  =  w,  /  ■JlRT  , 


and 


Q  =  PrS,  (4) 

in  which  N^2  (or  )  is  the  number  of  particles  to  be  transmitted  through  the  channel  from  the  inlet  (or  outlet)  to 
the  outlet  (or  inlet),  and  is  the  number  of  particles  coming  from  the  inlet  during  sampling  time.[3]  The 

quantity  Wj  is  the  x  component  of  the  bulk  velocity  crossing  the  open  boundary. 

EXPERIMENTAL  APPARATUS 

To  verify  the  present  numerical  results,  experiments  on  the  pumping  characteristics  of  a  HTDP  are  carried  out  by 
using  the  test  pump.  Dimensions  of  the  helical-grooved  rotor  of  the  test  pump  are  listed  in  Table  1.  The  rotor  is 
168.6  mm  in  diameter  and  132  mm  in  axial  length.  It  has  seven  threaded  rectangular  grooves. 

A  schematic  diagram  of  the  experimental  apparatus  is  shown  in  Fig.  4.  This  apparatus  mainly  consists  of  two 
parts:  the  test  pump,  and  instruments  that  facilitate  pressure  measurements  and  flow  supply  into  the  test  pump.  The 
test  pump  is  connected  to  a  two-stage  oil  rotary  pump  (970  l/min).  The  pressure  in  the  high-vacuum  side  is  measured 
with  a  Pirani  (7.6x  7.6  Torr)  and  Penning  (7.6x  1 7.6x  1 0"^  Torr)  gauge,  and  the  pressure  in  the  fore-vacuum 

side  is  measured  with  a  Pirani  gauge. 

Experiments  are  performed  by  varying  the  outlet  pressure  7^  in  the  range  of  0.02  Torr  <  <  4  Torr.  The 
rotational  speed  of  the  rotor  is  24,000  rpm  and  is  controlled  by  a  frequency  converter.  Test  gas  (iV^)  is  supplied 
through  a  mass  flow  controller  from  a  regulated  high-pressure  cylinder. 

RESULTS  AND  DISCUSSION 

The  present  numerical  results  are  compared  with  the  previously  known  numerical  and  experimental  ones 
obtained  by  Nanbu  etal.[2>],  as  shown  in  Fig.  5.  The  computations  are  conducted  by  employing  the  DSMC  and  N-S 
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methods  by  changing  the  inlet  pressure  for  the  fixed  outlet  pressure  =  0-3  Torn  The  present  DSMC  results  for 
Q3D  model  agree  well  with  those  of  Nanbu  et  al.  But  the  N-S  results  deviate  from  those  of  the  DSMC  method.  The 
throughput  Q  increases  linearly  as  the  pressure  difference  AP  {  =  P^-P^)  decreases.  The  DSMC  results  for  Q3D 
model  underestimate  by  as  much  as  20%  in  the  throughput  at  AP=  0.1  Torr  compared  to  the  experimental  data.  This 
discrepancy  is  partly  related  to  a  motion  between  the  rotor  and  casing  wall,  which  is  assumed  to  be  linear  in  Q3D 
model. 

Also,  the  DSMC  results  for  3D  model  are  shown  in  Fig.  5.  In  the  case  of  Q3D  model,  only  the  moving  effect  of 
the  casing  wall  on  the  performance  exists;  however,  in  the  3D  model,  Coriolis  and  centrifugal  forces  are  dominant 
effects  along  with  the  curvature  of  channels .  In  the  present  study,  these  effects  on  the  pumping  performance  are 
considered.  Consequently,  it  is  seen  that  the  DSMC  results  which  employ  3D  model  give  the  best  correspondence  to 
the  experimental  data. 

The  variations  of  /}  as  a  function  of  Q  at  ^2=  0.8  and  1.0  Torr  are  shown  in  Figs.  6  and  7  for  rotors  1  and  2, 
respectively.  The  value  of  P^  becomes  higher  as  Q  increases,  and  it  depends  on  .  Comparison  between  the 

experimental  data  and  the  DSMC  results  shows  good  agreement. 

The  effect  of  the  outlet  pressure  on  the  inlet  pressure  at  zero-flow  condition  is  shown  in  Fig.  8  for  rotor  1,  and  the 
comparison  between  the  experimental  data  and  the  numerical  results  is  also  made.  Comparison  between  the 
experimental  data  and  the  N-S  results  shows  good  agreement  in  the  range  of  P2>  2  Torr  (in  this  region, 
Kn  =  d ,  based  on  the  mean  free  path  at  the  outlet  and  the  channel  depth  «/  ,  is  in  the  order  of  10"^).  But  the 
discrepancy  between  the  two  results  becomes  larger  for  7^  -  114  Torr,  i.e.,  Kn  >  0.01 .  This  indicates  that  the  slip 

flow  analysis  of  the  N-S  method  is  inappropriate  for  Kn  >  0.01 . 

In  the  slip  flow  regime,  it  has  generally  been  recognized  that  the  layer  of  gas  adjacent  to  the  wall  is  no  longer  at 
rest  but  has  a  finite  slip  velocity.  Thus  the  continuum  assumption  is  not  valid  near  the  wall.  From  the  momentum 
transfer  between  gas  and  wall,  the  slip  on  the  wall  surfaces  is  given  by: 


+- 


2-0 


on 


Jw 


(5) 


where  a  is  the  accommodation  coefficient,  and  refer  to  the  slip  velocity  and  the  wall  velocity,  respectively. 


FIGURE  4.  Schematic  diagram  of  experimental 
apparatus. 


FIGURE  5.  Performance  curves  at  P2=03  Torr. 
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In  the  above,  n  is  the  non-dimensional  coordinate  perpendicular  to  the  channel  wall  that  is  normalized  by  the 
channel  depth  d  ,  and  subscript  W  indicates  that  the  variable  is  evaluated  at  the  wall. 

The  boundary  condition  of  Eq.  (5)  is  viewed  as  an  expansion  of  the  velocity  gradient  at  the  wall  in  powers  of 
Kn .  For  continuum  flow,  Kn-^0,  Eq.  §)  reduces  to  the  familiar  no-slip  condition.  Specifically,  the  no-slip 
condition  and  Eq.  (5)  are  the  zeroth  and  the  first-order  approximations,  respectively.  For  the  second-order 
modification,  the  following  slip  velocity  can  be  utilized  in  the  case  of  cr  =  1 : 


f 

Mj  +Kn 

V 


dn^ 


Jw 


(6) 


As  seen  in  Fig.  8,  the  results  obtained  by  the  N-S  method  with  second-order  SBC  are  nearly  the  same  as  those 
with  first-order  SBC  in  the  range  of  >  1  Torr.  The  differences  between  the  two  results  appear  noticeably  for  < 
0.1  Torr.  Although  the  N-S  method  with  second-order  SBC  fails  to  predict  quantitatively  the  experimental  data,  it 
seems  that  the  results  obtained  by  second-order  SBC  agree  qualitatively  with  the  experimental  data  compared  to 
those  obtained  by  first-order  SBC. 
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FIGURE  6.  Inlet  pressure  vs.  throughput  of  rotor  1. 
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FIGURE  8 .  Inlet  pressure  v^.  outlet  pressure. 


FIGURE  7.  Inlet  pressure  vs.  throughput  of  rotor  2. 
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FIGURE  9.  Performance  comparison  of  rotor  1  with 
rotor  2. 
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Experimental  results  are  also  compared  with  those  obtained  by  the  DSMC  simulations.  The  good  agreement 
between  the  two  is  shown  in  Fig.  8.  This  clearly  reveals  the  ability  of  the  DSMC  method  to  predict  the  molecular 
transition  flows.  Especially,  the  3D  model  predicts  well  the  experimental  data  compared  to  the  Q3D  model. 

In  general,  the  computing  cost  of  the  DSMC  simulation  is  linearly  proportional  to  the  gas  density  and  the  size  of 
a  computational  cell,  which  should  be  scaled  with  the  mean  free  path.  Therefore,  a  high  density  flow  requires  a  large 
number  of  computational  cells.  This  means  that  more  particles  and  more  collisions  between  molecules  must  be 
accounted  for  in  the  computation.  Consequently,  for  10"^  <  Kn  <  1,  although  the  DSMC  method  provides  more 
accurate  solutions  than  does  the  N-S  method,  it  requires  much  more  computational  effort  and  higher  expenses  as  Kn 
decreases. 

The  effect  of  the  outlet  pressure  on  the  inlet  pressure  at  zero-flow  condition  for  rotor  2  is  shown  in  Fig.  9.  For 
comparison,  the  experimental  data  for  rotor  1  in  Fig.  8  are  also  included  in  this  figure.  The  N-S  results  presented 
here  are  those  obtained  by  the  second-order  method.  It  is  seen  that  the  maximum  attainable  inlet  pressure  is  about 
6x10'^  Torr  for  rotor  1  and  4.3x10'^  Torr  for  rotor  2,  respectively.  As  indicated  in  Fig.  8  for  rotor  1,  the  N-S  results 
for  rotor  2  also  deviate  from  the  experimental  data,  except  the  range  of  ^  2  Torr.  On  the  other  hand,  as  expected, 
the  DSMC  results  give  good  correspondence  to  the  experimental  data. 

CONCLUDING  REMARKS 

Numerical  and  experimental  studies  have  been  made  for  the  molecular  transition  and  slip  flows  in  the  rotating 
helical  channel  of  a  HTDP.  The  computations  are  conducted  by  employing  the  N-S  and  DSMC  methods. 
Experiments  were  performed  by  varying  the  outlet  pressure  in  the  range  of  0.02  Torr  <  P^<  4  Torr. 

A  new  3D  model  of  a  helical  channel  has  been  developed,  and  the  effects  of  the  Coriolis  and  centrifugal  forces 
on  the  pumping  performance  are  discussed.  The  numerical  results  indicate  that  the  present  3D  model  gives  more 
accurate  solutions  than  does  the  previous  Q3D  model. 

Also,  the  present  results  obtained  by  the  N-S  and  DSMC  methods  agree  well  with  the  experimental  data  for 
.^K<0.01 .  In  particular,  the  N-S  method  with  SBC  drastically  simplifies  the  simulation  and  saves  computational 
effort  in  predicting  the  performance  for  this  regime. 
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